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Abstract 



ILJ We derive an analytic formula at three loops for the cusp anomalous dimension FcuspC^) 

*0 in AA = 4 super Yang-Mills. This is done by exploiting the relation of the latter to the Regge 

;-( limit of massive amplitudes. We comment on the corresponding three loops quark anti-quark 

potential. Our result also determines a considerable part of the three-loop cusp anomalous 
dimension in QCD. Finally, we consider a limit in which only ladder diagrams contribute to 
physical observables. In that limit, a precise agreement with strong coupling is observed. 
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1 Introduction 

The cusp anomalous dimension Fcusp (0) was originally introduced in [1] as the ultraviolet 
(UV) divergence of a Wilson loop with a cusp with Euclidean angle 0. It describes a wide 
range of interesting physical situations. One that will be of special importance in this paper 
is its relation to the infrared (IR) divergences of massive scattering amplitudes and form 
factors, see [2, 3] and references therein. 

In planar Af = 4 SYM, using dual conformal symmetry one can argue that the (Euclidean) 
Regge limit ^ s » m'^ , t of Coulomb branch amplitudes [4] is also governed by this function 



^We call this limit "Euclidean" because s and t have the same sign, as opposed to the usual Regge limit. 



Tcusp [5] • In this paper, we use this relation to extract the three- loop value of r^usp from the 
scattering amplitudes. 

The two-loop contribution to the cusp anomalous dimension in QCD was computed in 
ref. [6] and rederived and simplified in ref. [7]. In super symmetric theories, it is natural 
to define a generalization of the usual Wilson loop. As we discuss below, one can define 
a locally supersymmetric Wilson loop that couples to scalars in addition to gluons. The 
two-loop result for the latter loop in AA = 4 SYM is closely related to the result for the 
standard QCD loop, and was obtained in [8,9]. We reproduced these formulas as a check of 
our calculation, see eqs. (26) and (27). 

We find the previously unknown three-loop contribution to be given by 
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- 6 Li5(e-2'^) - 4 H2,3ie-'n " 6 H^,2{e-^^) - 6 H,^i(e-^^) + 4(3 U^ie'^'') 
-2C2Li3(e-2'^) + 2C2C3 + 3C5|. (1) 



Here ^ = tanh(y9/2, and Li„ are poly logarithms. Here </? = zc/) is the lorentzian version of the 
angle (a boost angle). The Hij are harmonic polylogarithms [10], whose definition we recall 
in section 2. This is one of the main results of this paper. 

There are several interesting limits of this function. 

• In the large angle cp ^ 00 limit, it grows linearly with the angle, FcuspCv') °^ V^r??isp- The 
coefficient T'^ is the anomalous dimension of a light-like, or null, Wilson loop. It is also 
related to the high-spin limit of anomalous dimensions of composite operators [11-13]. 
It is determined exactly by an integral equation [14]. 



• 



Small angle limit. At = 0, we have the straight line, which is 1/2 BPS, and its loop 
corrections vanish. The order 0^ correction to the BPS configuration is related to the 
energy loss of an accelerated quark. We have FcuspC^) = -(j)^ B{\,N) + 0{(f)^), where 
the "Bremstrahlung function" i?(A,A^) is exactly known [15,16]. In [15,16] this function 
was also related to other observables, such as the power radiated by a moving quark. 



the two point function of the displacement operator on the Wilson loop and the stress 
tensor expectation value in the presence of a Wilson line. 

• ^cuspi4') gives the quark anti-quark potential on S^ for a configuration which is sepa- 
rated by an angle 5 = vr - 0, see figure 1(b). The limit 6-^0 gives the quark anti-quark 
potential potential in flat space. That limit is considered in section 4. 

In M = A super Yang Mills we can also introduce a second angle at a cusp [17]. This 
second angle is related to the fact that the locally supersymmetric Wilson loop observable 
contains a coupling to a scalar. This coupling selects a direction n, where n is a point on S^. 
The Wilson loop operator is given by [18] 

W ~ TrlPe'^f A-dx+f \dx\n4i /2\ 

where we wrote it in Euclidean signature. One can consider a loop with a constant direction 
n, with drU = 0. Consider such straight Wilson line making a sudden turn by an angle 0, see 
figure 1(a). At the cusp we could consider the possibility of changing the direction n by an 
angle 6, cos6 = n-n', where n and n' are the directions before and after the cusp. In that 
case the Wilson loop develops a logarithmic divergence of the form 



-rcusp(<?i,6')iog 



uv 



{W)-e''''''^'''"''"^^^ (3) 

where Ajr/uv are the IR and UV cutoff energies, respectively. Thus, we have a cusp anoma- 
lous dimension, Tcuspi4',(^), which is a function of two angles and 6. The former is the 
obvious geometric angle and the latter is an internal angle. 

The same generalized cusp anomalous dimension Tcuspi4',(^) also characterizes the planar 
IR divergences that arise when scattering massive W bosons on the Coulomb branch of 
Af = 4 SYM. There, cos 6' = ni ■ n2 is the angle between the Coulomb branch expectation 
values ($) = dia.g(hini,h2n2,. ■ ■ ) associated to a pair of color adjacent external W bosons 
Wi^i and Wi^2- This generalized cusp anomalous dimension was computed to leading and 
subleading order in weak and strong coupling in refs. [17] and [9], respectively. We expect 
that, at three loops, Tcnspi'^,^) can be obtained from the result at 6* = of eq. (1) by setting 
^ = (cosh (f - cos 6)1 sinh (f. 

Ladder diagrams give the "most complicated" contribution at three loops. That is the ^^ 
piece in eq. (1). More generally, we find that up to three loops ladder diagrams are the only 
contribution to the ^^ term where L is the number of loops. We expect that to be true to 
all loop orders. One can therefore wonder if there is a physical quantity that is computed 
by summing only ladder diagrams. We find that this is indeed the case! Having a second 
angle 6, there are new limits one can take. In particular one can consider the limit where 
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(a) 

Figure 1: (a) A Wilson line that makes a turn by an angle in Euclidean space. The 
vectors n, n' denote the direction that sets the coupling to the scalar, (b) Under the plane 
to cylinder map, the same line is mapped to a quark anti-quark configuration. The quark 
and antiquark are sitting at two points on S^ at a relative angle oi 6 = tt - (p. Of course, they 
are extended along the time direction. 



i^ ^ oo, A -^ with A = Ae*^/4 (and therefore also X^^^) held fixed. That limit select the 
scalar ladders diagrams and is the subject of section 3. The sum of these ladder diagrams 
can be performed by finding the ground state of a one dimensional Schrodinger problem. It 
is possible to take the large A limit. We compare this with the strong coupling answer in the 
large e*^ limit and we find that they agree. In principle, it did not have to agree, since the 
order of limits is different. However, as in the BMN [19] or large charge limits [20], the two 
limits commute and we get a precise agreement. For a particular angle, namely = 0, we 
can solve it completely as a function of the coupling. 

There are two appendices. In appendix A we comment on direct relations between in- 
tegrals appearing in massive form factors and the Wilson loop under consideration. We 
explicitly give a sum of two three-loop Wilson loop diagrams that contributes directly to the 
three-loop form factor in QCD. We also compute a certain type of Wilson loop integral with 
interaction vertices (or equivalently, form factor integral) analytically to all loop orders. We 
observe that up to a simple factor, the answer at L loops is a homogeneous polynomial of 
degree 2L - 1 in the cusp angle (p and vr. Finally, appendix B contains the result for the 
Regge limit of individual four-point integrals that contribute to Fcusp- 
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Figure 2: Sample two-loop diagram contributing to the four-particle amplitude. Solid and 
wavy lines denote massive and (almost) massless particles, respectively. The precise masses 
are given by the labels. Dual conformal symmetry implies that the same function A4{u,v) 
describes two different physical situations: The Regge limit s ^ cx3 of (a) is equivalent to the 
Bhabha-type scattering (b), where the outer wavy lines have a small mass that regulates the 
soft divergences. 



2 rcusp(A,(/?) from scattering amplitudes 



2.1 Regge limit and soft divergences in massive AA = 4 super Yang- 
Mills 

The cusp anomalous dimension FcuspCv^) appears in (supersymmetric) Wilson loops with 
(Minkowskian) cusp angle ip. It can also be obtained from a leading infrared (soft) divergence 
of amplitudes or form factors involving massive particles [2]. 

In particular, it can be extracted from massive amplitudes on the Coulomb branch of 
AA = 4 SYM [5]. These amplitudes are obtained by giving a vacuum expectation value to 
some of the scalars of AA = 4 SYM. The string theory dual of this setup [21] suggests that 
the amplitudes defined in this way have an exact dual conformal symmetry [4]. Consider 
the four-scalar amplitude A^ with on-shell conditions p"^ = -(hi - /ij+i)^, where /ij+4 = hi. 
Here hi are four nonzero eigenvalues of the Higgs fields, which we take to all point in the 
same direction. The Mandelstam variables are s = (pi +^2)^ and t = {p2 +^3)^. ^ A priori 
the amplitude could depend on five dimensionless ratios built from the Poincare invariants 



^We follow the - + ++ metric conventions of ref. [4], so that s is negative for positive center of mass energy. 
The amplitude will be real for s and t both positive. 



s,t, hf. Dual conformal symmetry implies that it is a function of two variables only, 

M{s,t,hi,h2,h^,K)=M{u,v), (4) 

where 

_ hih3 _ h2hA . . 

u and V are invariant under dual conformal transformations. M.(^u^v) is known explicitly at 
one loop, where it is given in terms of logarithms and dilogarithms. It is an interesting open 
question what its functional dependence is at higher loops. In this paper, we determine the 
M ^ limit up to three loops. 

Following ref . [5] , there are two mass configurations that are of special interest two us: 

(a) The equal mass configuration hi = m, which implies that the external states are mass- 
less, p"^ = 0, and we have u = rn^js , v = rn^jt. In the limit m, v « 1, with u/v fixed one 
obtains the (IR-divergent) massless amplitudes at the origin of the Coulomb branch. 
In particular, in the four-point case discussed here, the latter are known to all loop 
orders. Here m acts as an IR regulator. 

(b) The two-mass configuration hi = h^ = m, h2 = h^ = M, where the on-shell conditions 
are pj = -(M - m)'^, and u = ni^js , v = M'^/t. In the limit of M » m, the kinematical 
configuration is that of Bhabha scattering. Here we have external massive particles of 
mass M , and m acts as an IR cutoff if we keep s and t of order M'^. In string theory, 
the different masses come from strings ending on stacks of D3 branes at different radial 
coordinates in AdS space. See ref. [4] for more details. In the planar limit we can view 
this as the scattering of a massive quarks and a massive antiquark. 

Thanks to dual conformal symmetry, both (a) and (b) are described by the same function 
M{u,v). 

We will be interested in the limit u « 1 of A4{u,v) [5]. This limit has different physical 
interpretations for cases (a) and (b). The two configurations (in the limit) are illustrated in 
Fig. 2. In the interpretation (a) of M{u,v), this is the Regge limit s ^ oo. In (b), in the 
limit we have heavy particles of mass M -m f^ M that interact by exchanging light particles 
of small mass m. The mass of the light particles regulates soft divergences. This is similar in 
spirit to giving the photon a small mass in QED. In the planar limit, and from the point of 
view of the low energy degrees of freedom, we are scattering a massive quark and anti-quark 
pair. Thus there are only two regions that give rise to soft divergences-^. See Fig. 2(b). The 



■^If we were scattering massive adjoint particles there would be four regions that give rise to soft diver- 
gences. 
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overall soft IR divergence can be obtained in the eikonal approximation. This leads naturally 
to a cusped Wilson loop [22]. One can then identify the coefficient of the IR divergence of 
the amplitude with (minus) the coefficient of the UV divergence of the cusped Wilson loop. 
Therefore we expect [5] ^ 

\og{M{u,v)) -^ (logM)reusp(A,</^) + 0(MO). (6) 

Note that there are two soft regions of the amplitude that contribute. These regions corre- 
spond to either the upper or lower propagator of mass m in Fig. 2(b) having soft momentum. 

We have the perturbative expansion 

r,,3p(A,y.)=^(-^) rSCv'), >^ = 9'ymN. (7) 

L>1 \o'^ / 

The logarithmic contribution on the r.h.s. of (6) comes from soft exchanges between the 
heavy particles. In this soft region, the amplitude can be approximated by form factors. In 
this limit, the Minkowskian cusp angle is a natural variable. It is defined as 

cosh 09 = — = , (8) 

where p2 and p^ are the (ingoing) momenta forming the cusp. The relation to the Euclidean 
cusp angle is = iip. Using the definition v = M^/(p2 + Ps)^ and the on-shell conditions 
P2=P3 = -M^, we have 

coshyj = 1 + — . (9) 

2v 

Another natural variable, which is also used frequently in the literature e.g. on Bhabha 
scattering integrals, is 

., -^-' -e-- (10) 

For < f < oo (i.e. in the Euclidean region) x is real and ranges from to 1. The inverse 
relation to v is 

X 1 

""^(T^^sinh^f ■ ^^^^ 

Notice that this relation is invariant under x -> l/x. ^ This symmetry corresponds to 
ip -^ -ip. In fact, we will find that our results, which are functions of v, are invariant under 
this transformation.*' 



Note that we changed the sign and normaUzation in the definition of rcusp(A, (p) w.r.t. ref. [5]. We have 

r*ii^r(^'V') = 72rX'(A,<p). 

^Note that in the physical region, where v has a small imaginary part coming from the Feynnian propagator 
prescription, the sign of this imaginary part has to be included in the x ->• l/x transformation. 

^ The fact that Tcusp should be analytic around a; = 1 or </) = is very clear when we interpret it as a 
quark- antiquark potential on the sphere, see figure 1. 



11 





13a 




Figure 3: Integrals contributing to the four-particle amplitude to three loops. Solid and 
wavy lines denote massive and massless propagators, respectively. Overall normalizations of 
s and t, as well as a loop-momentum dependent numerator factor in I^f, are not displayed. 



2.2 Evaluation of the four-particle amplitude to three loops 

The four-particle amplitude on the Coulomb branch has the perturbative expansion 

A 






(12) 



To three loops, it is given by [5,23,24] 



M'^^^ {U, V)=-- [haiu, V) + Isaiv, u) + 2hi,{u, v) + 2hi,{v , u)] . 



(13) 
(14) 
(15) 



Here /i, /2, /3a, -^36 are massive scalar four-point integrals. They are shown in Fig 3, and their 
precise definition is given in ref. [5], where they have been evaluated in various limits. Note 
that the expressions in eq. (13)-(15) are valid for generic M,f, i.e. for finite values of the 

mass(es) [24]. 

We are interested in the Regge limit n « 1 . Note that the rate of divergence of the scalar 
integrals without numerators can be predicted a priori. It follows from their topology and 
can be determined by the counting procedure of ref. [25]. For example, the one- and two-loop 
integrals have the following small u expansions 



h{u,v) — > \oguIi^i{v) + hfliv) + 0{u) , 

hiu.v) > \ogul2,l{v) + hfliyV) + 0{u) , 

h{v,u) -^ \og^ul^2ii^) + hgull^{v) + r2^o{v) + 0{u). 

u->-0 ' ' ' 



(16) 
(17) 
(18) 



(Here the superscript r indicates that the expansion terms in the last hne come from I2 
rotated by 90 degrees, i.e. with u and v exchanged.) For the integral Isb, which has a loop- 
momentum dependent numerator, the counting rules do not apply. Its Regge behavior was 
analyzed in ref. [5]. 

We wish to determine the coefficients of the logarithms in the Regge limit. In order 
to compute Fcusp at three loops, we need all lower-loop coefficient functions (including the 
0(u^)) ones, and the logu terms at three loops. Note that all 0{u) terms can be safely 
dropped. This is an advantage over dimensional regularization. 

The four-particle integrals in the Regge limit are very closely related to form factor 
integrals. One could in principle use the method of differential equations [26] in order to 
compute the latter. Here, we use a shortcut by making a convenient ansatz for the result. 
In order to do that, it is helpful to have an idea of the kind of functions that can appear in 
general. The solutions to the type of differential equations mentioned above can typically 
be written in terms of harmonic polylogarithms (HPLs) [10] of argument x. These functions 
are defined recursively by integrating the following kernels, 

fi{x) = -^, /o(x) = -, /_i(x) = -^. (19) 

I- X X 1 + X 

The starting point for the recursion are the weight one functions 

Hi{x) = -log{l-x), Ho{x) = log{x) , H_i{x) = log{l + x) . (20) 

Higher weight HPLs are defined in the following way, 

Ha,,a„...,aA^) = f fa^ Ha,,...,aM dt ■ (21) 



The subscript of H is called the weight vector. We use an abbreviation common in the 
literature. If m zeros are to the left of ±1, they are removed and ±1 replaced by ±{m + 1). 
For example, 

i^o,o,i,o,-i(^) = ^3,-2(x). (22) 

In related studies of massive vertex-type integrals, this set of functions was sufficient to 
describe all occurring integrals, see e.g. [27-29]. 

Motivated by this fact, we wrote down an ansatz in terms of harmonic polylogarithms 
for each integral, and determined the coefficients in this ansatz by evaluating various limits, 
for X -^ 0, or a; -^ 1, and by comparing against generic data points obtained from numerical 
integration. In doing so, we found the Mathematica implementation of harmonic polyloga- 
rithms [30] and the Mellin-Barnes tools [31] useful. Following this procedure, we were able to 
find analytical formulas for all coefficient functions needed to compute Fcusp at three loops. 

10 



As an example, at one loop, we have. 



/i,i =^21ogx, 

/i,o =e H' -4//i,o(x) -4if_i,o(a;)] . 



(23) 
(24) 



Here we introduced the useful abbreviation ^ = (1 - a;)/(l + x). The formulas are valid for 
< X < 1 and can be analytically continued to other regions, as we describe below. 

The results for the two- and three-loop integrals are given in Appendix B. We now proceed 
to present the result for Fcusp- 



2.3 Analytic three-loop result for Fcusp 

The perturbative expansion of the cusp anomalous dimension is 

A 



rcusp(A, ^) = E ( ^^ J riusUv') > > ^ = ^YM^ , 

L>1 VO^ / 



(25) 



where qym is the Yang-Mills coupling, and A^ the number of colors. As was already men- 
tioned, the result up to two-loop results was known. 

As explained above, we compute Fcusp (v^) by evaluating the (Euclidean) Regge limit of 
the four-particle amplitude, thanks to eq. (6). The Regge limit of all integrals contributing 
to eq. (12) can be found in Appendix B. 

We find the following results to three loops, valid in the Euclidean region a; > 0, 
1 



^(1) 

cusp 
^(2) 



(^) 



e[2iogx]. 



nuipi^) = -2^ 



rSsp(<^) = "2^ 



-- logo; (log X + vr^) 

^ log^ X + 2 logx (C2 + Li2(x2)) - 2Li3(x2) + 2(3] , 



logx (log^a; + n'^y 



- \ e[K2H^2fl{x) - AC3Ho,o{x) - 8C2H2,o{x) + 16C2i^-i,o,o(a;) - 40(2^^0,0,0(0;) 

- 16(2^^1,0,0(2;) + 32/7_4,o(a;) - 32H^^o{x) + 24iJ_2,o,o,o(2;) - 24iJ2,o,o,o(a;) 
+ 32/f_i,o,o,o,o(a;) - 88/7o,o,o,o,o(a;) - 32i/i,o,o,o,o(a;) - 25C4^o(a;) - I2C5' 

- ^ e[lQC3Ho,o{x) - 32(2/^-2,0(2;) + 32(2/^2,0(3;) + 64C2//o,o,o(a;) - l28H_,^o{x) 

+ 128i/4,o(a;) + 64i/_3,-i,o(a;) - 64i7_3,o,o(a;) - 64iy_3,i,o(a;) + 64iy_2,-2,o(a;) 

- 64i7_2,2,o(2;) - 64/72,-2,0(2;) + 64iJ2,2,o(a;) - 64/73,-1,0(2;) + 64/73,o,o(a;) 



(26) 



(27) 
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Figure 4: Analytic structure of T^nsp- The Euclidean region is x > 0. Below threshold, x is 
a phase, x = e^'^. Above threshold, we have that -1 < x < 0, with x having an infinitesimal 
positive imaginary part. The zigzag line denotes branch cuts along [-oo,0]. 



+ 64/73,i,o(^) - 64iJ_2,o,o,o(x) + G^H^aM^) + l92Ho,o,oflfi(^) - 32C3i/-2(x) 
+ 32C3i^2(a;) - 32C2//-3(a;) + 24(4^^0(0;) + 82(2^3(2;) + 8C2C3 + I2C5 



(28) 

8,9], 



The one- and two-loop results were known and confirm those quoted in the literature 
see also [6,7]. The three-loop result is new. 

Using simple relations between harmonic polylogarithms of argument x^ and x, and 
their product algebra (see e.g. [30]), we can rewrite the ^^ and ^^ terms at three loops in 
a simpler way. Moreover, replacing the HPLs Hn{x'^) appearing in the ^^ term by classical 
polylogarithms Li„(x^), as at two loops, leads to eq. (1) given in the introduction. This is 
the main result of this section. 

Let us discuss the branch cut structure of Fcusp- The formulas presented above are valid 
in the Euclidean region, < a; < 1. They are manifestly real for < x < 1. Some of the 
functions have branch cuts along x e [-00, 0] or x e [1, cx)]. Recall that we expect the results 
to have the symmetry x -^ 1/x. Therefore the latter branch cuts should be spurious. One 
may check this in the case of polylogarithms using the identity 

n 



Im(Li„(x + iO)) 



log"'' X ., 



X > 1 



(29) 



{n-iy. 

where iO indicates that x has an infinitesimally small positive imaginary part. In general, 
one can easily verify this property, as well as the x -^ 1/x symmetry, using relations between 
HPLs of argument x and 1/x [10,30]. 

In the physical region, we have to remember the iO prescription of the Feynman prop- 
agators. The sign of iO depends on whether we discuss the Wilson loop, or the scattering 
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amplitudes (and it depends on the metric conventions.) We will take x to have a small pos- 
itive imaginary part. There is a threshold which corresponds to the creation of two massive 
on-shell particles at x = -1 (or, equivalently, v = -1/4, i.e. t = -AM'^), which naturally divides 
the physical region into two domains, above threshold and below threshold. Below threshold, 
X lies on the unit circle in the complex plane, x = exp(z0), whereas above threshold x is real 
and satisfies -1 < x < 0. The analytic continuation to these physical regions is discussed 
in [28]. We summarize the analytic structure for complex x in Fig. 4. 

We stress that the final result for Fcusp is much simpler than the individual contributions 
presented in Appendix B. We already mentioned that there are different ways of writing the 
final result. More interesting than what specific functions are used to express the result in 
might be the question whether its symbol [32,33] has any special features. In fact our ansatz 
that the result could be expressed just in terms of harmonic polylogarithms implies that 
only the entries l±x, x can appear in the symbol. Inspecting the above results, we see that 
the symbol can be written in a way such that its entries are either x or (1 - x^)/x. ^ This 
property is refiected in the fact that we only needed harmonic polylogarithms with indices 
and 1 (when using x^ as argument). It is also remarkable that the index 1 appears at most 
twice. Finally, we remark that the first entry of the symbol is always x, as required by the 
branch cut structure of the loop integrals. 

Let us now consider various interesting limits of F^usp- 

• In the large ip limit (small v limit), Fcusp(A,v?) is linear in ip, and the coefficient of ip 
is the anomalous dimension of a light-like Wilson loop [6] . We reproduce the result for 
illustration, and as a check of our conventions'^, 

hrn F,,,p(A,y.) = Jy.F-,p(A) + 0{^') , (30) 

where 

rr..„(A) = 2 ( A) _ 2c. ( A)^ , no {^f . o(xn . (31) 

These results are correctly reproduced by eqs. (26) - (28). (In order to take the limit 
on the three-loop result, the simpler form given in eq. (1) is convenient, since its 
logarithmic dependence as x ^ is already manifest.) Moreover, since our result is 
exact in v, sub-leading terms in this limit can also be obtained. In fact, based on 



^The fact that the symbol can be written just in terms of the entries x and (1 - x'^)/x impUes that it has 
a symmetry under x -> -x. Note however that this is not a symmetry of the function. In particular, the 
latter acquires an imaginary part for x < 0, while it is real for a; > 0. 

®See footnote 4. We introduced a factor of 1/2 on the RHS of eq. (30) in order to have the same definition 
for r~gp(A) as in ref. [5], and elsewhere in the literature. 
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the discussion in [5], we might expect that a stronger version of eq. (30) should hold, 
namely 



1 



hm r,,,p(^) = - ^r-.p(A) + ^o(A) + 0{e-^) , (32) 



where 



is the collinear anomalous dimension in Higgs regularization. Equation (32) can be 
justified if one can show that the limits limt,«i lim„«i and lim^«i liim.u,v«i ,u/v Hxed com- 
mute for M{u,v). It was observed in ref. [5] that to three loops, these limits do not 
commute in general for the individual loop integrals contributing to M{u,v), however 
they do commute for M(u,v) to that order. 

Assuming that eq. (32) holds true, it is interesting that it relates the collinear anoma- 
lous dimension Qo{a) for massive amphtudes to the cusp anomalous dimension. 

In the opposite limit, small (p (large v), we find 

rSp(^)=f C4^' + o(^'). (34) 

Taking into account eq. (9) we see that in the limit l/v f^ ip'^. Hence we have, for the 
leading terms as y? ^ 0, 



rcusp(A,v5) = f 



2 4 / \ x3 



5(^)4(s^)^S(s^)-<^') 



+ 0{p'). (35) 



This expansion is in perfect agreement with the three loop expansion of the exact result 
in [15,16]. 

• The limit x ^ -1. The point x = 1 corresponds to the threshold v = -4M^ of creating 
two massive particles. It corresponds to 6 = tt -(f) ^ 0, which gives the quark antiquark 
potential in fiat space, see figure 1. In other words, when the quark and antiquark 
are very close to each other on the sphere, their potential is the same as the one they 
would have in fiat space. 

In order to reach this limit, we first have to analytically continue to the physical region, 
where x < (and has a small imaginary part coming from the Feynman prescription). 
We observe that the terms multiplying ^,^^,^^ in Fcusp vanish as x ^ -1 (reached after 
suitable analytic continuation, as discussed above). 
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For one and two loops, we verified the limit given in [9]. It is, in our notation (recall 
that X = e^^ ,(j) = iip, 6 = n - (p), 

di(0-7r) =-^ + ..., (36) 





ra(0--)=-|iog(^) + .... (37) 

At three loops, we find 



l-(3) /J, \ 8 4 1 1 

cuspVV^ J s d^ 5 



167rlog'(^ — j + yvr^logj^ — j 



/2d\ 4:7v^ 
+ 167rlog I — ) + 367rC3 + 247r 



(38) 



Most of these contributions come from the ^^ terms, though the ^^ term contributes 
with a I and -^ terms that have been already included in (38). The term of order 
C,^ does not contribute in this limit. The logs, as well as the 1/5^ seem to be in 
contradiction with the naive expectation. We will discuss the interpretation of this 
result in section 4. 



3 A limit that selects the ladder diagrams 

As we observed in explicit perturbative computations, the term with the highest power of 
^ seems to be given by ladder diagrams up to three loops. One can wonder if this pattern 
continues to all loops. Here we argue that it does and that there is a particular limit of 
rcMsp(0, 6*) that isolates such terms. 

Notice that such terms contain the highest powers of cos^. This highest power of cos^ 
can only come from a diagram where there are L scalars ending on each Wilson lines at the 
Lth loop order. The only such diagrams are ladder diagrams. We conclude that if we take 
the scaling limit 

A ^ , e^^ ^ oo , with A = — - = fixed (39) 

then the ladder diagrams are the only contribution to Tcusp — ^ ^iad(^, </>)• Note that after the 
scaling limit we get a non-trivial function of and A. In this section we derive a Schrodinger 
problem whose solution is the sum of ladder diagrams. We also solve it exactly for = 0, and 
for general angles at strong coupling. Notice that in this limit the coupling to the scalars 
takes the form Ze**^/^ + Ze^*^/^, where we have "+" for one of the lines and "-" for the other 
line. For the first line, in the scaling limit, we get 

V\- J dt{Ze''l^ + Ze-''l^)^Vl f dtZ (40) 
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For the other hne, the only surviving term is the couphng to Z. In this hmit all bulk 
interactions vanish and we have the free theory. So the problem is the same as having a free 
complex matrix field Z and computing the expectation value in the presence of the source 
(40), and the corresponding one for Z. 

3.1 Ladder diagrams and the Schrodinger problem 

It is convenient to think about the problem on the sphere. Then the scalar propagators are 

A C0S6' n/ \ -^ 1 /.IN 

^(^' ^) = ^T^? U ^ IT (41) 



Stt^ 2(cosh(r - 0") + COS0) ' Svr^ [cosh(r - a) + cos< 

where we wrote the propagator in the scaling limit (39). Here r and a are the global time 
coordinates of the two endpoints. 

As explained in [34] the sum of ladder diagrams can be performed by first introducing 
a quantity F(T, S) which is defined as the sum over all planar ladder diagrams, where we 
integrate all insertions up to r < T and a < S. Then the derivatives of F are 

drdsFiT, S) = P{T, S)F{T, S) (42) 

Then we write x = T - S and y = (T + S')/2 and we make an ansatz of the form F = 
^^g-Q„j/v]/^(x). Then eq. (42) becomes 



-S'^- ' ' 

'-^'T' 



Stt^ (coshx + COS0) 



ij(x) = -— ^(x) = 2EschH^) (43) 



In principle we could find all the eigenvalues of this problem. Note however, that we are only 
interested in the long time behavior of the sum, which is governed by the lowest eigenvalue, 
Qq. We then get Z ~ e^oC^™'^). The sign of Qq is not obviously fixed by (43), but we have 
fixed it here, with Qq > 0, so that get the expected sign at weak coupling, for example. 
Note that the Schrodinger problem, (43), has a discrete set of bound states with negative 
Schrodinger energy^, Esch, (real Q) and a continuum with Esch > 0, or fi imaginary. 

We remark that for given A and (p, it is possible to compute the ground state energy 
numerically. 

3.2 Exact solution for = 

In the particular case of = 0, the potential in eq. (43) becomes the Poschl- Teller potential 
cosh" (x/2) [35]. In this case the Schrodinger problem can be solved exactly by a variety of 



^ Since the potential is negative, there is always at least one bound state. 
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techniques. By the change of variables 



1 



1 + e^ 



one can map equation (43) to a hypergeometric equation, 



-1 + 



+ 1 



7r^ 



(44) 



(45) 



Imposing the decaying boundary conditions at x = ±oo we see that the possible eigenvalues 
for the bound states are 

Qn = ^o-n , 0<n<QQ (46) 

where n is an integer. This gives a finite number of bound states, which depends on the 
coupling. This is in contrast to the sum over ladders for the anti-parallel lines in flat space 
where the number of bound states is infinite beyond a certain coupling [36]. In particular, 
the ground state wave function is particularly simple 

1 



-ipo 



(47) 



(coshx/2)^o 

We can also easily compute the first correction away from the = limit by expanding 
the potential in (43). The first term is proportional to 0^/cosh {x/2). We can sandwich it 
between the ground state wave function (47) to get the first correction 



2flQ6fiQ 



A 



,(^o| 



(2cosh^x/2)' 



m 



A 



no(fio + l) 



4 47r2" (^oll^o) 47r2" 4(4^0(^0 + 2) + 3) 

Here we used the following formula in order to compute the expectation value, 

b 



(^oIlV^o) /(-I) ' 

Then the sum over ladders, up to order 0^ is 
1 - VkTI 02 



1 



(4J 



(49) 



■^lad 



1 + VT 



+ K 



+ 0(0') 



K, 



A 



vr^ 



(50) 



2 16 \1 + k + 2Vk + 1, 

One could also compute the 0^ term by using second order perturbation theory. 

We can use eq. (50) as a consistency check of our perturbative calculation of the ladder 
diagrams. The sum of the ladders at L loops is given by the ^^ term of Tculp that was 
computed in section 2. Expanding those terms for small angle, we find 



-^lad 



K 



4 24 



+ 0(0^) 



+ K^ 



16 288 ^ ^ 



+ K-' 



1 4302 
' 32 ~ 3456 



+ 0(4>') 



+ 0{k^) 



(51) 



This is in perfect agreement with the small k expansion of eq. (50). 
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3.3 Perturbative solution in A 



Here we show how to solve the Schrodinger equation for any angle perturbatively in the 
couphng A. 

To leading order at weak coupling we can approximate the potential by a delta function, 
since the energy is very small, 

1 2^ 



-5{x) 



(52) 



cosh X + cosh y? sinh ip 

which makes sure that we get the right result at first order in the coupling. 

To obtain the solution at higher orders, we use the following procedure. It is convenient 
to perform the change of variables 



\I'(x) = rj{x)e 



-nox/2 



(53) 



At leading order, we have rj^x) = 1, in agreement with eq. (52). 

Since the exponential factor is the correct solution as x -> oo, we can set the boundary 
condition //(oo) = 1. This normalizes the solution, but also is stating that we do not have 
the growing solution so that we have picked a unique solution of the equations. So, when we 
integrate the equation all the way to a; = we will find an ^o dependent value for the first 
derivative at the origin. We note that we can determine Qq from rj thanks to the boundary 
condition 



2dx\ogr](x)\:r=o = ^c 



(54) 



9.^(x)U = o - 

which follows from the x ^ -x symmetry of the problem and that the ground state wave- 
function is symmetric. Defining a new variable w = e'^ we have 



d^wdyjir) = -Vtody^Tj + k 



1 



1 



7] 



X 



Stt^ sinh (p 



_w + e'-P w + e '-P 

In this form, it is clear that the differential equation has four regular singular points, w 
0, oo,-e'^, -e^"^, so that it is a particular instance of the Heun equation. Expanding Vt^ 
kVLq + k^VLq + . . . and r] = 1 + kr]^^^ + . . .. At order A, eq. (55) becomes 



(55) 



This equation is easily integrated 



w + e'P w + e~^ 



(56) 



w w' 



1 



w" + ef w" + e'"^ 



(57) 



where we used the boundary conditions ai w = (or x = oo). In fact, in order to determine 
i^Q , we only need to carry out the first integration. We find 

n'^^^ = -2d^r]^'\^i=2^ (58) 

Remembering that the contribution to the ladders is given by F'"'^ = -Qq, we find perfect 
agreement with eq. (26), in the limit i6 — > oo (39). 

One can see that at any loop order L, the result for fig can be written as a (2L- l)-fold 
iterated integral, multiplied by a factor of (A/sinhy^)^. 

3.4 Comparison with strong coupling 

At strong coupling the potential becomes very deep and we can approximate the energy by 
simply the minimum of the potential, at x = 0. This then gives 

r^«'^ = -(]o = ^— — , for A»l (59) 

27rcos I 

We have verified that this is in agreement with the strong coupling computation of ref . [9] . 
Namely, [9] find that the cusp anomalous dimension at strong coupling has the form 

T,^,p-VXF{(P,e) (60) 

If we now take iO 

(61) 

We see that inserting (61) into (60) we get (59). In order to see this from the formulas given 
in [9], one sets q = ir with real r, and expands for small p. In this way one can check that 
the coefficient matches precisely with (59). 

This matching is a bit surprising. The ladders limit, A -> with A fixed, is different from 
the strong coupling limit A -^ oo with i6 » 1 fixed. That is, the result could in principle 
depend on the order in which the limits iO -^ oo and A ^ oo are taken. A heuristic explanation 
for this match is the following. As we discussed around (40), in the large e*^ limit one of the 
Wilson lines is sourcing mainly the Z field and the other the Z field. Thus the configuration 
can be viewed as many Z fields with a low density of derivatives. These derivatives build up 
the Wilson loop. 

The situation is therefore similar to the BMN limit [19], or more precisely the classical 
large charge limit discussed in [20]. In that limit a similar match with strong coupling is 
found for the first two orders in the effective coupling. Here we have checked the leading 

19 



oo. 


we 


find that 














F(0,^) 


ex 


gJ6l/2 




cosf 



order term. In principle, the formulas in [9] should also allow us to compute the subleading 
term. From the Schrodinger problem point of view, it is trivial to compute the subleading 
term, one simply has to consider the harmonic oscillator approximation around the minimum 
of the potential to get 



r''"^ = -^o = - '^^^ ^"^^ +0(1/71), for A»l (62) 

27r cos I 

and further orders could be computed by straightforward perturbation theory. 



4 The anti parallel lines limit 

When 6 = 7T - (j) ^ we have a quark and anti-quark on the sphere separated by a very 
small distance. In this limit we expect to find the quark-antiquark potential in fiat space. 
However, the approach to the limit is tricky because it does not commute with perturbation 
theory. Namely, if we first expand for small 5 for fixed A we get 

Tcusp{<l>,e,X) - -^S^i^ + 0{5') (63) 

where agat is the coefficient for a quark antiquark potential in flat space for two straight 
lines. That is, 

^fiat(r,A,g) = - """^^'^^ =-lim ;^log(H^D) (64) 

where W^ is a rectangular Wilson loop of separation r along a time T. On the other hand 
if we first expand in A and then, for each fixed order in A, we go to small 6, we do not find 
the naively expected behavior (63). The reason is that there is a mixing between the color 
degrees of freedom of the quark-antiquark pair and low energy degrees of freedom in the bulk, 
namely bulk modes of energy of order A/r. This causes IR divergences in perturbation theory 
starting at two loops in Af = 4 SYM [34,37] and at four loops in QCD [38].^*^ The origin 
of this is the following. A quark-antiquark pair in the singlet color combination can emit a 
massless adjoint and then be left in the adjoint. The force between the quarks in the singlet 
and adjoint channel are different. In the planar approximation, the force in the adjoint is 
zero. For example, a long time exchange of a scalar would contribute to the potential Vflat 
an expression of the form^^ 

f \ CO 

■ ^ ^ = ^/^(e"*^='"^'^'-l) (65) 
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'^'^Note that in the QCD hterature on the quark antiquark potential, the loop order is sometimes defined 
differently. 

^^The -1 subtraction in (65) is a counter term in the effective IR theory canceling a power UV divergence. 
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where l^ingiet = ■^ + 0{X^) is minus the singlet potential. This is the binding energy we loose 
when we emit an adjoint, in the planar approximation. Here 1/t^ is the large t approximation 
to the propagator. The UV cutoff of the effective IR theory, euv, is of the order of the 
separation r. 

This integral is perfectly finite in the IR. However, if we first expand in A, then we see 
that the term of order A^ has a logarithmic IR divergence, the term of order A^ has a power 
law IR divergence, etc. If we are on the sphere, we should replace (65) by 

( a b = — f (e-*^-ns_i) (66) 

27r2 J 2(cosht-l)^ ^ ^ ^ 

where now V^ingiet = 4^ + O(A^) and euv is of order 6. If we first expand the integrand in 
powers of A we now get integrals that are perfectly IR convergent. The long time cutoff is 
just the size of the sphere that we have set to 1. To first and second orders, namely, at two 
and three loops the diagram (66) contributes 



(B^r[-¥--] - (b^) 



- - ' ''''ry-\ 



to Fcusp, in agreement with the leading terms in (37) and (38). Here we have set euv ~ 5 in 
(66). Note, in particular, that the funny 1/5^ term is coming from a power law IR divergence 
in the three loop computation of the quark anti-quark potential in flat space. The log (5) 
term in the three loop result (38) can be obtained from the exponentiation of the two loop 
single log term in (67), see [37]. If instead of expanding in A we do the integral in (66) and 
then expand in 6 for S « X, then we get the expected behavior 

at two loops. In order to get the subleading coefficients represented by the dots one needs 
to be more careful about the matching of the effective theory in the IR with the full theory, 
see [37] for a full discussion^^. In the next section we will reproduce the two loop result 
of [37] by resumming ladder diagrams. We also do a similar resummation for the three loop 
ladder result. However, in the full theory the non-ladder diagrams do contribute. We leave 
the problem of performing the full three loop resummation for the future. 



^^The integral (68) has a logarithmic UV divergence, which is absorbed as part of the matching procedure 

[37]. 
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4.1 Ladders diagrams and the quark anti-quark potential 

Only ladder diagrams contributed to the one and two loop results (36), (37) of section 2. 
This means that a resummation of ladder diagrams would correctly reproduce the two loop 
quark anti-quark potential. Moreover, in section 3 we considered the limit iO -^ oo where only 
ladder diagrams survive. To compute the potential between the corresponding special quarks 
one has to re-sum these diagrams. Hence, in this section we will perform that resummation 
up to three loops. 

The Bethe-Salpeter equation maps the sum of ladder diagrams to the ground state energy 
of the Schrodinger problem [34] 

A 



" 47r2(x2 + l) 



(]2 



(69) 



as 



Qo = -r V{X) = Aci + A^cs + A^cg + . . . (70) 

The computation of Qq also contains the IR effects we discussed above. Therefore, the 
perturbative expansion in A for the Schrodinger problem (69) is a bit subtle. For example, if 
we applied the same iterative strategy we took in section 3.3 when solving for the Schrodinger 
problem on the sphere (43), we would get divergent integrals. Instead, we take a different 
strategy as we now explain. 

We divide the x axis into the large x (IR) region and the small x (UV) region. The two 
regions overlap at x ~ 1/V A where the potential and the energy are of the same order. In 
the IR region the equation can be solved straightforwardly. There, the solution up to two 
loops is a Bessel function. At three loops we also have to expand the potential to next order 
in the l/x expansion. In the UV region, we start from a constant wave function and iterate 
the equation to the desired order. The energy is then obtained by demanding that the two 
solutions match in the overlap region x ~ 

1 



A. In that way we find 



Cl 



C2 



C3 



1 

87r3 _ 
1 


log + 7e 1 
[flog ^.7.^ 


l^r.„ ^ ^^, ^A 


7 
~2~ 


7r2 


327r5 


[i^°^2vr'^^^jr^2vr^ '^ ' ^ 


' 12 



(71) 



where je is the Euler-Mascheroni constant. The two loop result C2 as well as the leading log 
in cs agrees with the computation of [37] using the effective field theory techniques discussed 
above. The two loop quark anti-quark potential at finite 6 is given by simply substituting 
A -^ A(l + cos^)/2. At three loops, for the first time, also non ladder diagrams contribute 
and the corresponding result differs from C3. We leave that computation to future work. 



22 



One can also study the Schrodinger problem in the limit where A ^ oo. In that limit we 
get a huge potential well and to leading order the ground state energy is simply given by the 
value of the potential at x = 0. That is 

no = — + 0{l) (72) 

n 

The exact same value is obtained by taking the iO ^ oo limit of the strong coupling result 
of [18]. That is a special case of the match obtained in section 3.4, in the limit -^ vr. 
In [34] the sum over ladder diagrams was compared to the strong computing computations 
in [18,39] for zero 6 and it was found to disagree. However, if one takes a quark antiquark 
with a relative 9 angle, which was also computed in [18], and one takes the large iO limit, 
then one finds that it matches precisely with (72). 

At small A, the potential in (69) has a single bound state. As we increase A more and 
more bound states go down from the continuum. Beyond a critical value A > Ac = vr^ we 
have an infinite number of bound states [36] . ^^ One may wonder whether also these can be 
matched with the modes of the string at strong coupling. We find that unlikely as generically 
the string has infinitely many degrees of freedom and not just one. Note however that the 
density of states near the top of the potential at Vtoo - does match with the string density of 
states for the relevant mode [36]. This can be interpreted as the statement that the classical 
motion for the string mode identified in [36] matches the classical (or large A) limit of the 
Schrodinger problem. 

5 Conclusions 

In this paper we have performed a three loop computation of the cusp anomalous dimension 
for the locally supersymmetric Wilson loop in AA = 4 super Yang Mills. Equivalently, we have 
computed the three loop quark/anti-quark potential on the three sphere. The final result 
is eq. (1). The computation was done by considering the scattering of massive W bosons 
and focusing on the IR divergences, which are given by rcMsp(v?). In turn, this amplitude is 
related by dual conformal symmetry to the Regge limit of the four point amplitude of massless 
particles in the Higgs regularization [4,5]. We have discussed some checks of this result. In 
particular, we have matched the leading 0^ term to the exact computation in [15, 16]. We 
have also discussed the limit that corresponds to the quark anti-quark potential. In this 
limit, small 5, we get a result that is different from the 1/5 behavior expected naively, where 



^•^In [36] it was also argued that this should be the behavior in the full theory for the antiparallel lines in 
flat space. 
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S is the relative angle. This is due to some divergences that appear in perturbation theory 
in this computation [34,37,38]. After taking into account the long distance effective theory 
describing these IR effects one can explain the behavior of the result for small 6. In principle, 
one should be able to do a resummation and give the full three loop result for the flat space 
quark-antiquark potential. This is a problem we left to the future. 

We have also identified a large i6 limit (39) of the generalized cusp anomalous dimension, 
rcMsp(0, 6'), where only ladder diagrams contribute. These ladder diagrams can be summed 
by solving a simple Schrodinger problem. For general angles the potential does not appear 
to be solvable, but one can develop a simple perturbative expansion where one clearly sees 
that the answer is given by iterated integrals that give Goncharov polylogarithms. For 
a particular angle, = 0, we could solve the problem and sum the ladder diagrams for all 
effective couplings A. It is also possible to get the first term in the 0^ expansion around = 0. 
The final answer is in (50). Also, we have taken the small S limit of the Schrodinger problem 
and expanded the answer to third order in A. We reproduced the two loop resummed answer 
in [37]. It is very simple to compute the strong coupling limit of the result for the ladder 
diagrams. More precisely, we take the effective coupling A to be very large. Interestingly we 
find a precise match with strong coupling string theory computations. Here one first takes 
A to be large, and then takes e*^ -> oo. In this limit the classical string theory answers go 
like e*^/^ so that multiplied by vA one gets V A which is the behavior of the strong coupling 
ladder diagrams. This agreement was not preordained, since the scaling limit (39) is explicitly 
pushing us towards the small A region. It seems that the origin of this agreement has the 
same underlying reason as the agreement for the BMN [19] or large charge limits [20]. 
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A Results for Wilson loop integrals relevant for QCD 

It is important to realize that, at weak coupling, Wilson loops depend only mildly on the 
particle content. There are two sources of differences between J\f = 4 SYM and QCD for 
these Wilson loops. The first is the additional coupling of the loop to scalars in AA = 4 
SYM. The second is the specific particle content, which e.g. enters at two loops through the 
one-loop gluon propagator correction. We want to stress that many diagrams that have to 
be calculated are identical. This means that our result computes a part of the QCD result. 

Moreover, there is a direct relation between some of the Wilson loop integrals we are in- 
terested in, and the infrared divergent part of form factor integrals discussed in the literature. 
These connections are very helpful. First of all, they give an idea about the kind of functions 
that can appear in such calculations. In some cases the relations are even more concrete, and 
they provide two different ways of computing the same object. We will give two examples. 
The first one relates the ^^ term in Fcusp to a known form factor integral. Inversely, we point 
out that our analytic result at three loops, obtained from scattering amplitudes, implies that 
some three-loop Wilson line integrals that are relevant for QCD are now known. In the 
second example, we discuss the horizontal ladder diagram, which contributes to the ^ term. 
This is just one contribution to the ^ term but there are other contributions which are not 
these ladder diagrams. Here, the Wilson loop picture turns out to be extremely simple, and 
allows us to derive an all-orders result for this class of integrals. Of course, to get the QCD 
answers we need to set ^ = ^r^. 

A.l Results for crossed ladder integrals at two and three loops 

It is easy to see from the Wilson loop computation that the term proportional to ^^ at two 
loops comes from the ladder with two rungs and a contribution from the squared one-loop 
ladder. It is well known that this combination results in a crossed ladder Wilson line integral, 
see Fig. 5(a). The relevant UV divergent part of that integral can be equivalently obtained 
from the IR pole of the form factor integral of the same topology. The latter is given by 
V6/4mi;-i in Tcf. [27]. Comparing to eq. (27), we find perfect agreement. 

Having obtained the full three-loop result for Fcusp by evaluating the four-point amplitude 
in the Regge limit as described in section 2, we can provide an analytical answer for some of 
the Wilson loop diagrams, and equivalently, form factor integrals that are of more general 
interest. In fact, the ^^ terms in Fcusp come from summing ladder diagrams only. Taking into 
account the exponentiation of lower-loop graphs, they correspond to the sum of the Wilson 
loop diagrams shown in Fig. 5(b,c). Therefore, we have that 

[Fig. 5(6) + Fig. 5(c)] - log A x \^^ term in eq. (1) + C(A°) , (73) 
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(a) (b) (c) 

Figure 5: Crossed ladder diagrams that appear in the computation of Fcusp at two and three 
loops. 



where A is some UV cutoff. Similarly, in dimensional regularization log A would be replaced 
by 1/e. Being built from pure gluon propagators, these diagrams also appear e.g. in the QCD 
calculation of Fcusp at three loops, or equivalently when computing massive form factors at 
three loops. 

A. 2 Exact result for a class of diagrams 

Here we compute the Regge limit of the horizontal ladder diagrams, see Fig. 6(a), to all loop 
orders, by relating them to known off-shell ladder diagrams [40]. We do not know if there 
is a physical limit that isolates only these diagrams. For massive scalar integrals, there is 
a useful counting procedure to determine their rate of divergence in the Regge limit [25]. 
Applied to the present case, we find the following behavior, 

4oL„talladder(«>^) ^, ^OgU R^^\^) + 0{u') . (74) 

We want to compute the function R^^\ip) that multiplies the logarithm. It is well known 
that there are simplifications for that term at the level of Feynman parameter integrals [25] . 

Here, we follow a different approach. We find it more convenient to think of the limit 
M ^ as the soft limit in the Bhabha scattering picture, as discussed above. At L loops, it 
is clear that the leading soft divergence comes from regions where all loop momenta are soft. 
There is a factor of 2 because there are two relevant regions. In this limit, it is easy to see 
that the ladder becomes a three-point diagram with eikonal lines, cf. Fig. 6(b). The latter 
can be represented as a double line integral, cf. Fig. 6(c), and we obtain a result proportional 
to (for L > 1), 
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(a) 



(b) 



(c) 



Figure 6: Horizontal ladder diagram (a) in the soft limit m^/s -^ 0. Thick solid lines denote 
heavy particles of mass M, thin solid lines denote particles of mass m, and wavy lines denote 
massless particles. The leading IR divergence of (a) is given by twice or IR divergence of 
the vertex integral (b). The divergent part of the latter is captured by the eikonal limit, 
leading to the Wilson line integral (c). The doubled lines in the latter denote the Wilson 
line contour. At L > 1 loops, the massless coordinate space integral in the interior can be 
identified with the off-shell ladder integral $(^^1). 



where 



y'lin) 



nP2: 



1/2(^-2) = r2P^, yi2 = yi-y2, 



(76) 



and where $(^) are known off-shell ladder integrals. They are given by a linear combination 
of logarithms and polylogarithms of homogeneous degree 2L [40] . 

By taking the eikonal limit, we have in fact introduced UV divergences for small ri 2, and 
as a result, the integral in eq. (75) is formally zero. The coefficient of the IR logarithm we 
want to compute is given by minus the coefficient of the UV logarithm. The latter can be 
easily extracted from the small ri 2 region of the integral, by introducing a cut-off for the 
radial integration. Collecting all proportionality factors, we find 



i?(^)(¥,) 



^ cosh ip -1 f dz ^ 



rd^^(L-i)( 



+ 2z cosh ip + z'^' 1 + 2z cosh ip 



+ z^)' 



(77) 



where $ are rescaled off-shell ladder integrals, 

l>(^) (x, y) = A(x, i/)$(^) (x, y) , A(x, y) = y/{l - x - y^ - 4xy . (78) 

For L = 2, eq. (77) has appeared in the Wilson loop computation of ref. [8] and was evaluated 
analytically in ref. [9]. Here we show how to perform the integration for any L. 
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The integral can be easily evaluated by using the following integral representation for 
^^^^ [40], valid for x >y, 

Plugging this formula into eq. (77), we notice that the change of variables t = pC^ ,z = pe^ 
makes one integration trivial (after expanding powers of logarithms into sums), while the 
remaining one reduces to a known integraP"^. In this way, we arrive at the result 

(80) 

a;=0 



R^'^^'\v) = ig{L) 



/ d y''smh{x(p) 
\dx) sin(x7r) 



Here g{L) is given by a sum, 

2-2L+3 L-1 IL-1\IL-1\ '^L-l-rn2 

L!(L- 1)! .^^^2=oV mi )\ m2 )l + mi + m2 
and hence it evaluates to rational numbers at any L. Inspecting eq. (80), (81) we see that 

p(L.i). ^ , (-l)^^H47r)^^^\ (7r-z^\ 

where Bn{x) is the Bernoulli polynomial. Up to the ^ factor, the answer for R^^^ is a 
polynomial in vr and ifl The polynomial is homogeneous, and of degree 2L - 1. 

For example, at two and three loops, we immediately reproduce eqs. (85) and (94), 
obtained by an independent computation using Mellin-Barnes methods. At four loops, we 
obtain 

i?(4) ((^) = -^ J_ (31 ^ TT^ + 49 ip'^ TT^ + 2lip^n^ + 3 <^^) . (83) 

Let us perform a consistency check on this result by taking the limit x ^ 0. In the latter, 
^ ^ 1, and (f = -logx stays invariant. We can see that this is consistent with eq. (A.l) of 
ref. [41], where this limit was previously computed. 

At L loops, we know from unitarity cuts that these ladders appear with coefficient 
(-A/(167r2))^ in the amplitude. Therefore, they give the following contribution to T^usp 

^(honzontal ladders) ^ ^ ^_^j'' ^(i)^^^ 



i^ ^ \^ 



St^---.(^)- <-) 



An i^, (2L - 1) 



^''Gradshteyn and Ryzhik, "Table of Integrals, Series, and Products", Academic Press, Sixth Edition, 
p. 549, equation (14). 
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In summary, we have computed the leading Regge hmit of the horizontal ladder diagram 
shown in Fig. 6(a), or equivalently, the coefficient of the pole term in the corresponding 
massive form factor integral shown in Fig. 6(b), at any loop order L. Up to an overall factor 
of ^ = (cosh(y9 - l)/sinhy9, the result is a homogeneous polynomial of degree 2L - 1 m Lp and 

TT. 

Interestingly, the horizontal ladders result (82) also appears in the TBA integrability 
approach to the cusp anomalous dimension [42]. There, it appears as the leading Llischer 
contribution with L + 1 scalar insertions at the cusp. 

The computations are superficially different. Let us now give an argument for this agree- 
ment^^. The calculation of the Llischer contribution would be the following. Considers a 
Wilson loop with two cusps, and with Z^ and Z^ inserted at the two cusps, respectively. 
Here Z = $5+2$6 are complex scalar fields. The Wilson lines connecting the cusps are smooth 
and couple to scalars different from $5 and $6- At tree level, one only has Wick contractions 
between the Z and Z fields. We are interested in the ffist wrapping-type diagrams, where a 
field is emitted from one Wilson line, interacts with the Z-Z propagators, and is absorbed 
by the other Wilson line. This term will be linear in ^. We will focus on the term where the 
scalars are emitted and absorbed. We now argue that this calculation would lead to a result 
proportional to (82). 

The argument consists of two steps. First, notice that there is only one type of Feynman 
diagram we need to consider. The reason is that the S'f/(4) charges of the scalars have to 
flow through the diagram. The only way to connect the diagram in a planar way is to use 
four-scalar vertices. We are interested in the scaling dimension of the cusp with operator 
insertion, so we need to compute only the UV divergences of this diagram. The latter come 
from the integration regions close to the cusps. We can focus on one cusp to compute them. 
When doing so, at leading order the propagators connecting to the other cusp factor out, 
and we are left with the diagram of figure 6(c). This shows that the two calculations have 
to agree, up to the overall coefficient. 

One may wonder whether there is a situation where only the infinite class of integrals 
computed here contributes to a physical quantity. In section 3, we saw that this was the 
case for a similar class of ladder diagrams, which could be singled out by taking the large 
i9 limit, where 9 determines how the Wilson lines couple to the scalars. More generally, we 
can think of the scattering amplitude depending on several angles on S"^, as explained in the 
introduction. Then, the limit of section 3 just selects the vertical ladder diagrams. It seems 
possible that a similar limit selects the horizontal ladder diagrams. 
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They agree up to a minus sign. 
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B Results for individual integrals up to three loops 

Here we present analytic formulas for the integrals contributing to Fcusp to three loops. The 
Regge limit of the one-loop integral was already given in the main text. At two loops, we 
find 

I2,i=^y7r^hgx+-hg'x , (85) 

/2,o =e [ - 16//-3,o(2;) - lQH^,o{x) - WH2,o,o(x) - 8Ho,oM^) 

- 24C2Ho,o(x) - 4C,Ho{x) - 42^] , (86) 
Il,=e[2\og'x], (87) 
III =e' [ - lQC2Ho{x) - 16H2fl{x) - 16i/_i,o,o(a:) - 8i7o,o,o(x) - l6H,^o,o(x) - 4(3] , (88) 
^2,0 =e' [s2C2H^ifl{x) + 24C2i/o,o(a;) + 82(2^^1,0(2;) + l6H^3,o{x) + l6Hs,o{x) 

- 32H^2,-ifiix) + lQH^2fi,oix) + 32if_i,2,o(x) + 32/7i,2,o(x) + 16if2,o,o(a;) 

+ 32H2,i,oix) + 32/7_i,_i,o,o(x) + 16iJ_i,o,o,o(2;) + 32iJ_i,i,o,o(2;) + 8Ho,o,o,oix) 
+ 32iJi,_i,o,o(2;) + 16/f 1,0,0,0 (^) + 32if 1,1,0,0 (a;) 

- lQC2H.2{x) + SCsH.^ix) + 8C3ifi(x) + 32C2i/2(x) + 6OC4 •] • 

At three loops, we have 



haiu,v) 
ha{v,u) 



u^O 



u^O 



M^O 



logu/3„,i + C(nO), 

log^ U l3a,3 + log^ U J3"„,2 + log U Jg"^,! + C(m°) 
log^M/3f,,2 + logM/3fe,i + C(m°) , 



We found the following coefficient functions, 

ha,i = e ^ [Tvr^ log X + IOtt^ log^' x + 3 log' x] , 
which is in agreement with the calculation for general L, c.f. eq. (80), 

-^3a,3 =7.i.h,l) , 

Ila,2 =('[- 40C2i^o,o(x) + 16i7_3,o(x) - 48/f3,o(x) - 64/72,o,o(x) 

- 48i7_i,o,o,o(x) - 64ffo,o,o,o(a;) - 48i7i,o,o,o(a;) - 8C3i^o(^)] , 
Ila,i <' [32C3i/-i,o(a;) + IQGHo^oix) + 32C3i/i,o(x) - 32(2/^-2,0 (a;) + 224C2i/2,o(^) 



(89) 

(90) 
(91) 
(92) 
(93) 

(94) 

(95) 
(96) 
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+ 160C2i/-i,o,o(a;) + 256C2i^o,o,o(a;) + 160(2^^1,0,0(2:) + 256H^,o{x) - 64ff_3,-i,o(x) 
+ 128/7_3,o,o(x) - 64if_3,i,o(2;) - 6AH^2,-2,o{x) - 64/7_i,_3,o(a;) + 192i7_i,3,o(x) 

- 64i7i,_3,o(x) + 192/Ji,3,o(x) - 64/^2,-2,0 (x) + 256/^2,2,0 (x) - 64i/3,-i,o(x) 

+ 256/^3,0,0 (a;) + 192/73,1,0(3;) - 192i/_2,-i,o,o(x) + 160/f_2,o,o,o(a;) + 64/J_2,i,o,o(a;) 
+ 256/J_i,2,o,o(2;) + 256ifi,2,o,o(a;) + 64if2,-i,o,o(2;) + 2881/2,0,0,0(2;) + 320if2,i,o,o(a;) 
+ 192/7_i,_i,o,o,o(x) + 256i7_i,o,o,o,o(x) + 192i7_i,i,o,o,o(x) + 272iJo,o,o,o,o(x) 
+ 192ifi,_i,o,o,o(a;) + 256i7i,o,o,o,o(a;) + 192iJi,i,o,o,o(a;) - 16(3^-2(0;) + 48(3^^2(0;) 

- 96C2H^^ix) + 20AC,Ho{x) + 160(2^^3(2:) + 82(2(3 + 12(5] , (97) 

and for I^i,{u,v)\ 

hb,2 = + ^h,ih,i , (98) 

/3b,i =e [ - 128/^4,0(3;) - 64i/_3,o,o(a;) - 96ff3,o,o(a;) - 144^2,o,o,o(a;) - 128i7i,o,o,o,o(a;) 

- 216i7o,o,o,o,o(a;) - 32(2i72,o(a;) - 176(2i/o,o,o(a;) - 64(2iJi,o,o(x) 

- 16(3ifo,o(x) - 152UHo{x) - 24(5] . (99) 

Finally, we have 

^36,1 =e[^log'x + ^7r2log3x+^7rMogxl. (100) 



We remark that the factorization of the leading terms Jg 2, I^a 3' -^36,2 can be understood from 
the analysis of [25,41], where the systematics of the Regge limit were investigated. 
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